Quantized magnetic moment at the edge of a single-walled carbon nanotube 
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We investigate many-body effects near the edge of a single-walled carbon nanotube and find it 
turns magnetic with quantized edge moment solely depends on the chiral vector, i.e. the topology of 
the carbon nanotube. The distribution of the edge moment reveals an approximate supersymmetry 
even in a realistic model. Our findings demonstrate crucial importance of the interplay between the 
electronic correlation and the edge morphology. The edge moment provides an excellent candidate 
for nanomagnet farbrication which has potential applications in biological and chemical detections, 
and in spintronics. 



A single-walled carbon nanotube (SWNT) [1 can be 
viewed as rolling up the two-dimensional graphene [2H5] 
into cylindrical shape with radius of the order of nanome- 
ters. Its quasi-one-dimensionality leads to pronounced 
many-body effects, manifest in optical transitions for 
electron- hole excitons [6j [7] , electronic bandgaps by scan- 
ning tunnelling spectroscopy [8] etc. and leads to exotic 
ground states including Luttinger liquids [9] [10] , Wigner 
crystal [IT] and Mott-insulating spin liquids [12] . 

Because the open boundary in the honeycomb network 
changes the electronic states dramatically, it is interest- 
ing to explore how correlation effects reshape the physical 
properties. For the simplest zigzag edge, both analytic 
and numeric approaches [I3j [14] indicate additional zero- 
energy excitations at the edge, while for the armchair, 
none of these evanescent modes exist. However, count- 
ing the total number of the zero modes for generic chiral 
SWNTs remains unknown. Since these edge states are 
degenerate, inclusion of mutual interactions between elec- 
trons induce ferromagnetic coupling and promotes mag- 
netism in SWNTs. In fact, correlation-induced magnetic 
moment has be speculated for quite a while [15-18 but 
its magnitude hasn't been quantified yet. 

The central goal of this Letter is to show how edge 
moment emerges in generic SWNTs. Cutting a single- 
walled carbon nanotube perpendicularly, we investigate 
the many-body effects near the hydrogenated edge and 
find the nanotube turns magnetic with quantized edge 
moment solely depends on the chiral vector, i.e. the 
topology of the carbon nanotube. Meanwhile, the spin 
density for the edge moment only shows up in one of the 
sublattices for the honeycomb, revealing an approximate 
supersymmetry [TIB [20] even with the inclusion of the re- 
alistic band structures. Our findings demonstrate that 
the interplay between the electronic correlation and the 
edge morphology is of crucial importance at nanoscale. 
The edge moment, solely depending on the topology of 
the nanotube, provides an excellent candidate for nano- 
magnet farbrication which adds a spin flavor to biolog- 
ical [21] and chemical detections [22] and also potential 
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FIG. 1: A SWNT and its flattened form with chiral vec- 
tor (10,2), where the red and green spheres denote C and 
H atoms respectively. The spin density (grey color) arisen 
from electron-electron interactions appears near the edge and 
decays rapidly into the bulk. 



applications in spintronics [23H26] for carbon nanotubes. 

Let us start with how the topology of a SWNT is 
classified. The chiral vector, composed of two integers 
(A/ - , M), denotes a lattice vector in the honeycomb lattice 
of graphene and label the waist uniquely in the roll-up 
SWNT. To remove the redundancy from the hexagonal 
symmetry in the honeycomb structure, we fix the con- 
vention for the chiral vector, N > M > [27 . For 
M = 0, this family of SWNTs is named "zigzag" while 
the other known "armchair" family carries the chiral vec- 
tor, N = M. Otherwise, they are called chiral in general. 
To classify all possible edges of semi-infinite SWNTs is 
more tricky. For a given chiral vector, there are infinite 
ways to cut the nanotubes. Therefore, we only consider 
the primary edge by cutting along the chiral vector, i.e. 
the perpendicular cut around the nanotube. For simplic- 
ity, we assume the primary edge is terminated by hydro- 
gen atoms. As an example, a SWNT with chiral vector 
(10,2) and its flattened form are shown in Fig. 1. The 
honeycomb structure consists of two sublattices A and B. 
The primary edge is obtained by cutting the SWNT per- 
pendicularly. The unit cell contains 12 carbon atoms on 
each sublattice and the hydrogenated edge is connected 
to 2 carbon atoms on the A-sites and 10 on the 5-sites. 
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FIG. 2: Graphical solutions for z in the complex plane for 
SWNTs. Varying zigzag C — to armchair C = 1, each 
family of SWNTs falls onto a universal curve and forms the 
set of the eye-ball contours. As an example, the red dots 
represent the non-unitary solutions \z\ ^ 1 for semi-inifite 
(13,0) zigzag SWNT. 



We also show the spin density determined by our first- 
principles calculation, which will be explained later in 
detail; the spin density arisen from electron-electron in- 
teractions appears near the edge and decays rapidly into 
the bulk. It is surprising that the spin density profile has 
nodes on the A-sites, hinting an approximate supersym- 
metry in SWNT, even with the inclusion of realistic band 
structures. 

We start with the Hubbard model on the honeycomb 
lattice, 



h= J2 W4, 
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where the hopping amplitude t xy = —t for nearest neigh- 
bors and vanishes otherwise. The electron-electron inter- 
action is approximated by the on-site repulsion U > 0. 
Although this Hubbard Hamiltonian is a simplified ver- 
sion for the realistic SWNT, it is still too hard to be 
solved exactly. 

Ignoring the on-site repulsion momentarily, the best 
strategy to pinpoint the edge states in the presence 
of the open boundary is by the hidden supersymmetry 
(SUSY) pi [20] in the honeycomb lattice. The SUSY 
not only enforces these edge states at zero energy but 
also guarantees their wave functions vanish in one of the 
sublattices (A-sites in Fig. 1). The non- vanishing wave 
function on sublattice B can be solved by the generalized 
Bloch theorem 20 , 
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where ^ is a (TV + M) dimensional column vector de- 
scribing the all 5-sites in the unit cell, z is the eigenvalue 
of unit displacement operator and n labels the number 



of the unit cells in the longitudinal direction. The (— l) n 
factor is for notational convenience and q are appropriate 
constants to satisfy the boundary conditions at the edge. 
Since the translational invariance is lost in the longitu- 
dinal direction, in addition to the (\z\ = 1) plane- wave 
solutions, TV< evanescent modes with \z\ < 1 appear. As 
is clear in Fig. 1, there are M/N constraints from A/B- 
sites at the edge. However, the wave function has nodes 
on A-sites, rendering the TV constraints trivially satisfied. 
Subtracting the lost degrees of freedom from the M con- 
straints, we end up with TV e = TV< — M zero modes due 
to the edge of the SWNT. 

The solutions of Z{ for the evanescent modes can be 
found as in the standard calculations for the band struc- 
tures. The characteristic equation from the generalized 
Bloch theorem turns out to be remarkably simple, 



z-\ 



N 



1 



1 



M 



1. 



(3) 



Introducing the complex function f(z) = (z — 1) N (1/z — 
1) M , the above equation is equivalent to wcgf(z) = 
(modulus 2tt) and \f(z)\ = 1, where the latter constraint 
gives a set of "eye-ball" contours as shown in Fig. 2. 
Note that these contours only depend on the chirality 
C = M/TV, varying from zero (zigzag) to one (armchair). 
All contours intersect at two points at z± = e ±Z7r / 3 , which 
correspond to the pair of the Dirac conducting channels. 
All evanescent modes requires \z\ < 1 and thus must 
occur on the shorter segment between z± = e ±27r / 3 . Fur- 
thermore, one can show that axgf(z) changes monoton- 
ically along the eye-ball contour. As a result, counting 
edge states is equivalent to tracking the winding number 
of aigf(z) passing through zero when going from z+ to 
Z-. The factor [z — 1) N in the complex function gives 
phase change 2tt7V/3 while (l/z-l) M = (z-l) M (-z)~ M 
gives 2ttM/3 - (-2ttM/3) = 4ttM/3. Therefore, the to- 
tal phase change is A[aigf(z)] = 2tt(N + 2M)/3. For 
metallic SWNTs where TV — M is an integer multiple of 
three, straightforward analysis gives the winding number 
from the total phase change is 7V< = (TV + 2M)/3 — 1. 
Similarly, for semiconducting ones, the winding number 
is 7V< = [(TV + 2M + l)/3], where [x] is the Gauss sym- 
bol taking the greatest integer less than or equal to x. 
Subtracting the number of constraints M, 
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(semiconducting). (5) 



To investigate the correlation effects, let us consider 
the semiconducting SWNTs in weak coupling limit first. 
The effective theory in low-energy limit is described by 
the zero modes at the edge since the bulk is gapped. 
By the projection operator P e onto the edge states, 
the semi-infinite carbon nanotube is thus mapped into 



the finite A^-cluster problem. Since all edge states 
are pinned at zero energy, the hopping Hamiltonian af- 
ter projection vanishes. Only the interaction survives, 
H e = U ^2 X P e [n x ^n x ^]P e > 0. Consider the fully polar- 
ized state with all edge states in the spin-up states. It is 
clear that the fully polarized state is an eigenstate of H e 
with zero eigenvalue. Since the projected Hamiltonian 
is positive definite, it is the desired ground state with 
non- vanishing edge moment, 
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The weak-coupling result can be directly generalized to 
arbitrary interaction strength by the Lieb's theorem [28] 
which ensure the ground state is unique up to the trivial 
(25 + l)-fold degeneracy. 

The metallic SWNTs are complicated by the pair of 
one dimensional conducting channels from the Dirac 
cones. We can follow similar procedure to map out the 
gapped states in the bulk and derive the effective the- 
ory for the remaining edge states and the ID conduct- 
ing channels. Again, the on-site interaction after pro- 
jection prefers fully polarized edge states. However, due 
to the interaction between the edge states and the con- 
ducting channels, the magnitude of the total spin can not 
be determined. One may worry that Kondo-like screen- 
ing from the itinerant carriers can further suppress the 
edge moment, or even make it vanish. Evaluating the 
exchange integral between the edge and itinerant states, 
we found the exchange coupling is ferromagnetic, which 
becomes irrelevant under renormalization group transfor- 
mation [29]. Therefore, it is natural to expect that the 
edge moment will survive but its magnitude is beyond 
the reach of the weak-coupling approach. 

So far, we show that the quantized magnetic moment 
in a semiconducting SWNT only depends on its chiral 
vector. However, it is not clear whether this beautiful 
quantization is realistic or an artifact of the simple Hub- 
bard model. Therefore, we also carry out first-principles 
calculations within the local spin-density approximation, 
that takes account of the realistic band structure and 
the Couloumb interaction. The self-consistent band 
structure calculations under lattice optimization were 
performed using the full-potential projected augmented 
wave method [30] as implemented in the VASP pack- 
age [31 . As an illustrating example, we compute the 
spin density profile in the (10, 2) chiral SWNT as shown 
in Fig. 1. The semi-inifite nanotube was simulated by 
a finite segment, with one of the edges being the de- 
sired shape while with the other being the armchair (not 
shown in the figure) without additional edge states. We 
have checked carefully that our results do not depend on 
the length of the segment and no extra moment is spotted 
near the armchair edge either. When summing all spin 
densities on the B-sites, the total moment turns out to be 
M = 3/i#, as predicted in Eq. [6] Furthermore, the nodal 













FIG. 3: Edge moments in zigzag carbon nanotubes of various 
radii. The spin density distributions for (a) N = 5 (b)N — 6 
(c) N = 7 (d) N=8 (e) N =9 (f) N =10 are shown in grey 
color while the red and green spheres denote C and H atoms 
respectively. Except for the small nanotubes (a) and (b), the 
magnetic moment is localized near the edge. 



structure of the spin density is manifest and serves as a 
direct evidence for the approximate SUSY in nanotubes 
even with the inclusion of realistic band structures. 

We also explore how the moment changes as the ra- 
dius of the nanotube varies. Concentrating on the sim- 
plest zigzag SWNTs, the computed spin density profiles 
for various radii are shown in Fig. 3. For the smallest 
radius N = 5, no spin polarization is found. For N = 6, 
the moment starts to appear but is not localized to the 
edge. For larger radii where the curvature effect is neg- 
ligible, the edge moments are localized. As long as the 
moment is localized to the edge, the magnetic moment, as 
shown in Fig. 4, follows the quantization rule in Eq. ^ 
even for the metallic SWNTs! For semiconducting ones 
N = 7, 8, 10, the spin density profiles, shown in Fig. 3(c), 
(d), (f), decay rather quickly into the bulk and agrees 
with the predictions from the Hubbard model. For the 
metallic one N = 9 in Fig. 3(e), the spin density extends 
a bit further into the bulk due to the exchange coupling 
between the edge and the itinerant states but the mag- 
netic moment remains quantized. The Lieb's theorem 
seems to provides a natural explanation for the quan- 
tized moment within the Hubbard model approximation. 
But, it remains an open question why the quantization 
is so robust even in a more realistic model, as seen in the 
first-principles calculation. In addition, our numerical re- 
sults show that the hidden SUSY remains valid to high 
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FIG. 4: Magnitude of magnetic moment in zigzag SWNTs. 
The green filled diamonds represent the predictions from an- 
alytic calculations for semiconducting nanotubes. The red 
circles mark the numerical results from the first-principle cal- 
culations. Both approaches give robust quantization of the 
edge moment except for small circumferences in the shaded 
blue regime. For N < 6, the curvature effect can not be 
ignored and ruins the edge moment. 



precision as the spin density exhibits the nodal structure. 

In comparison with the nonmagnetic state, the energy 
gain of the spin quantized CNTs with N=7~13 range 
from 0.12 to 0.30 eV per cell, respectively, which cor- 
respond to about 20 meV per edge carbon for all the 
cases studied. To understand the deviations at N = 5, 6, 
we studied a corresponding fictitious system of a flat 
graphene ribbon with periodic boudary condition along 
the chiral vector, by the same numerical method. Al- 
though the system is topologically equivalent to the nan- 
otube, we find that the spin density is localized and the 
quantization of the total spin is restored by the fully spin 
polarized flat edge bands. Therefore, we conclude that 
the deviations found for nanotubes with smaller radii 
arise from curvature effects which also agree with pre- 
vious studies [27] on the band structure in the literature. 
It is straightforward to generalize the numerical calcula- 
tions for other chiral SWNTs as well. 

We have applied both analytic and numerical methods 
to establish the emergence of the quantized magnetic mo- 
ment near the edge of a SWNT. The edge moment can 
be detected (and perhaps manipulated) by spin-polarized 
scanning tunnelling microscopy [32] or by magnetic reso- 
nant force microscopy [33] . Because of the robustness of 
the edge moments in the SWNTs, it provides an excellent 
method to fabricate nanomagnets. Our findings demon- 
strate that the interplay between the electronic correla- 
tion and the edge morphology is of crucial importance 
at nanoscale. The edge moment, solely depending on 
the topology of the nanotube, provides an excellent can- 
didate for nanomagnet farbrication which adds a spin 
flavor to biological [21] and chemical detections [22] and 
also potential applications in spintronics [23-26 for car- 
bon nanotubes. On top of the surprising phenomena for 



correlated physics in one dimension, carbon nanotubes 
can be used for biological and chemical detections at 
nanoscale as well. The presence of the edge moment may 
enhance the performance with spin resolution. In recent 
experiment, it was demonstrated that the spin and the 
orbital motions of electrons in carbon nanotubes are cou- 
pled strongly. Thus, when electron passing through the 
edge moment, the spin-orbit coupling provides another 
knob to control the transport properties and may pro- 
vides another window for spintronics applications. 

We acknowledge Leon Balents and Matthew Fisher for 
valuable discussions. We acknowledge supports from the 
National Science Council in Taiwan through Grant No. 
NSC-97-2112-M-007-022-MY3. The hospitality of KITP 
in Santa Barbara is also greatly appreciated. 



[1] 

[2] 
[3] 
[4] 

[5] 

[e; 

[7] 
[8] 
[9 



[10 
[11 

[12 



[13 
[14 

[is; 
lie; 

[17 

[is; 

[19 
[20 
[21 

[22; 



J. C. Charlier, X. Blase and S. Roche, Rev. Mod. Phys. 

79, 677 (2007). 

K. S. Novoselov et a/., Science 306, 666 (2004). 

K. S. Novoselov et a/., Nature 438, 197 (2005). 

Y. Zhang, Y.-W. Tan, H. L. Stormer and P. Kim, Nature 

438, 201 (2005). 

A. K. Geim and K. S. Novosselov, Nature Mat. 6, 183 

(2007) and references therein. 

F. Wang, G. Dukovic, L. E. Brus and T. F. Heinz, Science 
308, 838 (2005). 

G. Dukovic et a/., Nano Lett. 5, 2314 (2005). 

H. Lin et a/., Nature Mat. doi:10.1038/nmat2624 (2010). 

M. Bockrath, D. H. Cobden, J. Lu, A. G. Rinzler, R. E. 

Smalley, L. Balents and P. L. McEuen, Nature 397, 598 

(1999). 

H. Ishii et a/., Nature 426, 540 (2003). 

V. V. Deshpande and M. Brockrath, Nature Phys. 4, 314 

(2008). 

V. V. Deshpande, B. Chandra, R. Caldwell, D. S. 

Novikov, J. Hone and M. Bockrath, Science 323, 106 

(2009). 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusak- 

abe, J. Phys. Soc. Jpn. 65, 1920 (1996). 

M. Fujita, K. Nakada, G. Dresselhaus and M. S. Dressel- 

haus, Phys. Rev. B 54, 17954 (1996). 

K. Wakabayashi, M. Fujita, H. Ajiki, and M. Sigrist, 

Phys. Rev. B 59, 8271 (1999). 

S. Okada and A. Oshiyama, Phys. Rev. Lett. 87, 146803 

(2001). 

T. Hikihara, X. Hu, H.-H. Lin and C.-Y. Mou, Phys. Rev. 

B 68, 035432 (2003). 

Y.-W. Son, M. L. Cohen and S. G. Louie, Nature 444, 

347 (2006). 

B.-L. Huang, S.-T. Wu and C.-Y. Mou, Phys. Rev. B 70, 

205408 (2004). 

T. Pereg-Barnea and H.-H. Lin, Europhys. Lett. 69, 791 

(2005). 

D. A. Heller, E. S. Jeng, T.-K. Yeung, B. M. Martinez, 
A. E. Moll, J. B. Gastala, and M. S. Strano, Science 311, 
508 (2006). 

E. S. Snow, F. K. Perkins, E. J. Houser, S. C. Badescu 
and T. L. Reinecke, Science 307, 1942 (2005). 



[23] T. Ando, J. Phys. Soc. Jpn. 69, 1757 (2000). 

[24] D. Huertas-Hernando, F. Guinea and A. Brataas, Phys. 

Rev. B 74, 155426 (2006). 
[25] K. C. Nowack, F. H. L. Koppens, Yu V. Nazarov and L. 

M. K. Vandersypen, Science 318, 1430 (2007). 
[26] F. Kuemmeth, S. Hani, D. C. Ralph and P. L. McEuen, 

Nature 452, 448 (2008). 
[27] N. Hamada, S. Sawada and A. Oshiyama, Phys. Rev. 

Lett. 68, 1579 (1992). 
[28] E. H. Lieb, Phys. Rev. Lett. 62, 1201 (1989). 



[29] A. C. Hewson, The Kondo Problem to Heavy Fermions, 

Chapter 3, Cambridge University Press (1997). 
[30] G. Kresse and J. Joubert, Phys. Rev. B 59, 1758 (1999). 
[31] G. Kresse and J. Furthmiiller, Phys. Rev. B 54, 11169 

(1996). 
[32] O. Pietzsch, A. Kubetzka, M. Bode and R. Wiesendan- 

ger, Science 292, 2053 (2001). 
[33] D. Rugar, R. Budakian, H. J. Mamin and B. W. Chui 

Nature 430, 329 (2004). 



